The solution of a chiral random matrix model with complex 

eigenvalues 



G. Akemann 



Service de Physique Theorique* CEA/Saclay 
F-91191 Gif-sur-Yvette Cedex, France 



Abstract 

We describe in detail the solution of the extension of the chiral Gaussian Unitary Ensemble 
(chGUE) into the complex plane. The correlation functions of the model are first calculated for a 
finite number of N complex eigenvalues, where we exploit the existence of orthogonal Laguerre poly- 
nomials in the complex plane. When taking the large- N limit we derive new correlation functions 
in the case of weak and strong non-Hcrmiticity, thus describing the transition from the chGUE 
to a generalized Ginibre ensemble. Applications to the Dirac operator eigenvalue spectrum in 
QCD with non-vanishing chemical potential are briefly discussed. This is an extended version of 



hep-th/0204068 



PACS 12.38.Lg, 11.30.Rd, 05.40.-a 



SPhT T02/053 



* Unite associee CNRS/SPM/URA 2306 



1 Introduction 



The application of Random Matrix Models (RMM) in a variety of physical systems has a long history 
0. Roughly speaking its applications can be divided into two classes. RMM can be either used to 
describe local properties. The correlations for example of energy levels often shown a universal behavior 
on the distance of the mean level spacing, making a RMM description possible. Choosing the RMM 
with the same global symmetries as the underlying Hamiltonian of the system fixes completely the 
form of the local spectral fluctuations. On the other hand RMM can be used as a combinatorial tool 
to count classes of graphs on random surfaces, which may then be interpreted as Feynman graphs of a 
corresponding Quantum Field Theory. Our investigations here belong to the first class of applications 
on a local level, and the corresponding correlation functions are often called microscopic. 

In most cases RMM have to be seen as effective models owing their validity from symmetry 
considerations alone. However, in some cases the mapping to the theory one wishes to describe can be 
made exact, as for example in dynamical triangulations of random surfaces and its relation to Quantum 
Gravity [0] or in the application of RMM to the Dirac operator spectrum in four dimensional QCD ||. 
At the present stage the RMM we propose here as a model for QCD with chemical potential is just 
an effective or phenomenological model. But a more precise mapping to QCD via chiral perturbation 
theory may also be possible. 

In many applications models of random matrices having real eigenvalues are considered. All such 
possible models have been classified j|] and all the corresponding correlation functions are known 
(for a complete list of results and references see e.g. H). On the other hand there are also physical 
situations where complex eigenvalues occur. As examples we mention localization in superconductors 
||, dissipation and scattering in Quantum Chaos or Quantum Chromodynamics (QCD) with a 
non- vanishing chemical potential ||] . Other possible applications include the fractional Quantum Hall 
Effect as in || or a two-dimensional plasma of charged particles JTC) , 1 1 1 . A classification of all models 
with non-Hermitian matrices has been achieved more recently [1^] and the number of possible cases 
exceeds by far those with real eigenvalues. But even when considering the simplest extension of the 
three RMM classes introduced by Dyson |l3|] much less is known. They consist of real symmetric 
{(3 = 1), complex Hermitian {(3 = 2) and quaternion self-dual matrices {(3 = 4) where the Dyson 
index (3 labels the ensemble by the power of the Jacobian diagonalizing the matrices, which is given 
by the Vandermonde determinant A. Due to their invariance these ensembles are called Gaussian 
Orthogonal, Unitary or Symplectic Ensemble, the GOE, GUE and GSE respectively. 

In a seminal work the extension of these RMM into the complex plane has been obtained by Ginibre 
|14|. He calculated the joint probability distribution function (j.p.d.f.) for (3 = 2 and 4 and partially 
for /3 = 1 as well as all microscopic correlation functions for (3 = 2. Since then progress has been slow. 
The full j.p.d.f. for (3=1 was calculated much more recently in but the correlation functions 
are still lacking up to date. The correlation functions for (3 = 4 were found in All these models 
are Gaussian ensembles with independently distributed matrix elements. More recently generalized 
Gaussian ensembles were considered containing a parameter r that permits to interpolate between 
real and fully complex eigenvalues. The global density interpolating between the Wigner semi-circle 
on the real line and the constant density on a circular disc of Ginibre was found in [17] (see also 
figl for numerical results of all three ensembles) . Together with the introduction of the interpolating 



parameter r the regime of so-called weak non-Hermiticity was discovered [|1JJ, 20] for (3 = 2 and 1. 
For = 2 it interpolates smoothly between the know microscopic correlations of the GUE and those 
of Ginibre called strong non-Hermitian limit |j~9| ]. The techniques to derive such results are either 
supersymmetry or orthogonal polynomials [21, 11] which also exist on the complex plane ||. 



For (3 = 4 the T-dependent correlations were discovered very recently in [22| in both the weak and 



strong limit. New correlation functions for (3 = 2 were also found in [23] where so-called Dirac mass 
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terms were introduced in the application to three-dimensional QCD with complex eigenvalues. 

From the RMM @ which was first proposed to describe the QCD Dirac spectrum of real eigen- 
values, the chiral Gaussian Unitary Ensemble (chGUE), we know that the model we are aiming for 
should belong to the class of chiral models. In the presence of a quark chemical potential it is also 
known that the Dirac eigenvalues become complex ||. We thus wish to define and solve such a chiral 
RMM with complex eigenvalues. The model is found by taking the know analogy between the GUE 
and the chGUE on the level of real eigenvalues and translating it into the complex plane. The starting 
point is the known complex GUE introduced in Q. In analogy to the real case and from the finding 
of Hermite polynomials in the complex plane for the non-chiral complex GUE in ||] we expect and 
do find Laguerre polynomials in the complex plane. This enables us to provide explicit solutions at a 
finite number N of eigenvalues as well as in the weak and strong non-Hermitian large- N limit. 

The situation we wish to describe, QCD in the presence of a chemical potential, is quite delicate. 
The difficulties encountered when trying to do first principle calculations by using Monte Carlo simu- 
lations on a lattice are reviewed in 25]. The are due to the problem of convergence for a determinant 
of a Dirac operator with complex eigenvalues. Despite the fact that recent progress was made in |26|| 
the full phase diagram with fi ^ remains unexplored so far. The analytical knowledge of microscopic 
correlation functions from complex chiral RMM could thus be very useful, in view of its success in 
predicting real eigenvalue correlations. 

The chemical potential problem has already been addressed in quite a few works using RMM. In 
H the failure of the quenched approximation was illuminated by Stephanov. He also calculated the 
global eigenvalue density and its envelope as a functions of \x. We will use his results to determine our 
non-Hermiticity parameter r as a function of fi for small values of /i. In [^Tj the temperature-density 
phase diagram of QCD with two light flavors was predicted using a RMM. A more refined model 
incorporating also vector and axial symmetries of the action was studied in [28]. Effective models for 
a phase transition with the spectral density vanishing to an arbitrary power were introduced in |3(J . 
Furthermore, the analyticity properties of the partition function [31], problems of convergence [32] as 
well as the global symmetry of such Dirac operators [^] were analyzed using RMM. Finally first data 
for quenched simulations of QCD with fi ^ have been confronted with RMM predictions, finding 
a crossover from the Wigner-Dyson class over strongly non-Hermitian, Ginibre type correlations to a 
Poisson distribution |33|]. There the nearest neighbor spacing of eigenvalues in the bulk from non-chiral 
RMM were compared to the data. Here, we wish to provide analytic information on all microscopic 
correlation functions at the origin, which are important for chiral symmetry breaking 34], using a 
chiral eigenvalue model. 

The paper is organized as follows. In section § we introduce a chiral complex eigenvalue model. A 
justification is given why to call it chiral and why in particular it represents a complex extension of 
the chGUE. The identification of our parameter r as a function of the chemical potential [i in QCD is 
made by comparing with the RMM pi . In the next section || we introduce the orthogonal polynomials 
belonging to our model, show that they are given by the Laguerre polynomials in the complex plane 
and give the exact solution for all correlators for a finite number of eigenvalues N. Sections || and [| 
treat the large- N limit at weak and strong non-Hermiticity, respectively. New correlation functions 
are found in both cases by taking appropriate limits of the Laguerre polynomials. We illustrate our 
findings with a few examples in each case. We also prove that the weak correlations interpolate 
between the known correlations of the chGUE and those in the strong limit from section [|. This is 
done by letting the weak non-Hermiticity parameter a approach zero or infinity in sections [| and [| 
respectively, explicitly checking that the correlations coincide. Our conclusions are presented in the 
last section || A brief summary of our results has been already published in ][35|] . 
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2 The model 



We first define our model in terms of complex eigenvalues Zj=i n and then explain in detail why it 
is the generalization of the chGUE into the complex plane. The partition function is defined as 



r N 



(2.1) 



where we have introduced the following weight function 



w 



(a) 



M 2a+1 exp 



N 



1 



f 2 i * 

2 {Z +Z 



(2.2) 



and An(zi, . . . , zn) = Y\k>i( z k ~ z i) is the Vandermonde determinant. The parameter a > — 1 can be 
taken real for the time being. In the application to QCD it will be an integer, a = Nf + is, where Nf 
labels the number of massless quark flavors and v is the sector of topological charge. The parameter 
t G [0, 1] controls the degree of non-Hermiticity. For r = 1 we are back to the chGUE with real 
eigenvalues, as the weight functions ( |2.2j ) then reduces to a 5-function in the imaginary part 9m(z) 
times a Gaussian weight for the real part $le(z). On the other hand for r = the degree of non- 
Hermiticity is maximal, and we will see the individual role of r in the weight of the real and imaginary 
part below. 

The motivation for the definition of the model ( |2.1| ) comes from a comparison to the non-chiral 
model of complex eigenvalues, the complex extension of the GUE, introduced in 0, 19, 23]. Let us 



therefore recall the relation between the GUE and chGUE for real eigenvalues first. The chGUE is 
defined in terms of a complex matrix W of size N x (JV + v) 



Z 



chGUE 



dWdWUet(WWY f e- N9tT(WWt) = f°° f[ (dA^f'+V 2 ^) \A N (X U . . . , \ N )\ 2 

Jo j= i 

/ II (dxjlx^f+^e-^) \A N (xl,...,x 2 N )\ 2 . (2.3) 

J — CO -_i 



In the first step we have diagonalized and given the partition function in terms of the positive eigen- 
valuesQ Xj of the matrix W while in the second step we have simply substituted Xj = x 2 -. Note the 
squares inside the Vandermonde due to this transformation. Now we compare to the GUE given in 
terms of a Hermitian matrix H of size N x N, 



Z, 



(2N f ) 
GUE 



dHdet(H) 2N fe 



-Ngtv(H^) 



N 



, ^ 2N 
(XX j Xj 



\&n(xi, 



,XN 



(2.4) 



Both models can be solved using the method of orthogonal polynomials |T(| and the corresponding 
sets are given by the generalized Laguerre polynomials for the chGUE and by the Hermite polynomials 
for the GUE (at Nf = 0, for Nf > see pq] ). Moreover, when taking the parameter a = Nf + v to be 
real we can reconstruct the orthogonal polynomials of the chGUE on the full real line from the even 
subset of the orthogonal polynomials of the GUE. Because of the squares in the Vandermonde in eq. 
( p.3j ) we only need even polynomials there. This relation also holds for more general, non-Gaussian 
weight functions (see e.g. in |p9| ). For a = —1/2 this is nothing else than the known relation between 
Hermite and Laguerre polynomials |37| (see eq. (|3.8|) below). 



1 To be precise these are called singular values. 
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Let us now turn to the complex generalization of the GUE as introduced in || IS, 23]. The 
Hermitian matrix H is replaced by a complex matrix 



J 



H + i 



1 + T 



-A 



(2.5) 



with t € [0, 1] . The weight for both the Hermitian and anti-Hermitian part of J, the Hermitian 
matrices H and A respectively, is chosen to be Gaussian with variance (1 + r)/2N. We therefore 



obtain |T1 



dJdJ^ det( JJ ] ) N f exp 



I Yl (dzjdz* \ Zj \ 2N fexp 
J 3=1 V 



1 — T 

N 



(2.6) 



1 -T 



2 V I ^3 



2 / 2 i *2- 

~ 2^ + *J < 



|Ajv(^i, ••• ,z N ) 



For averages over the matrix elements of J this corresponds to 



(JklJlk) 



T 

N 



{JklJkU 



1 

N 



(2.7) 



which illustrates the interpolating role of r: for r = 1 we have Jik = and thus Hermitian matrices 
whereas for r = all Jki are mutually independent corresponding to maximal non-Hermiticity. In 
the latter case the model was already defined and solved by Ginibre in [14] where the orthogonal 



polynomials are monomials in the complex eigenvalues Zj of J. Here, we can study the transition 
between the GUE and the Ginibre ensemble and the corresponding polynomials were found to be 
Hermite polynomials in the complex plane ||. Using the same reasoning as in the mapping between 
the real eigenvalues model of the GUE eq. ( |2.4| ) and the chGUE eq. ( |2.3| ) , eq. ( |2.6j ) immediately leads 
to our model eq. (|2.1| ) above. This way we have solved the problem of how to extend the positive 
real eigenvalues ( |2.3[ ) into the complex plane. We may of course ask ourselves to what kind of matrix 
representation eq. (2.1) would correspond to. We would have to express eq. (2.1) in terms of a 
different matrix J, probably with two different complex matrices now, in the form J ~ WW^ + iVVl 
However, we have not been able to find such a representation so far. 

Let us compare to a different matrix model which has been introduced by Stephanov for QCD 
with non- vanishing chemical potential Q, 



dWdW^ det 







iW + ii 




N r 



-NgtT{WW' r ) 



(2i 



Here, the matrix inside the determinant replaces the QCD Dirac operator with the chemical potential 
added as 70/x. The disadvantage from adding a constant matrix to the chiral random matrix is that 
the model can no longer be written in terms of eigenvalues of the matrix W. Our idea is to use the 
model eq. (2.1) instead and to suppose that on the scale of the mean level spacing of the complex 
Dirac operator eigenvalues the microscopic correlations are the same. A somewhat similar observation 



has been made in [38] where a chiral matrix model with temperature was considered. In order to 
capture the effect of the lowest Matsubara frequency a constant matrix i^ot was added to the Dirac 
matrix. An analysis in terms of eigenvalues of the matrix WW^ is also no longer possible, although 
the original Dirac eigenvalues remain real in this case. However, the microscopic correlations remain 
the same compared to the model at t = 0, when translating to the same mean level spacing f3£|. To 
avoid misunderstandings let us stress that we do not claim that the correlations at /j, = and ^ / 
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are the same but that (|2,8| ) at fi ^ and our model at r < 1 share the same microscopic correlation 
functions. 

In order to relate these two parameters we compare the macroscopic spectral densities of the two 
models. In our model it can be read off from jl7]] since the preexponential factors in the weight eq. 
( [2.2| ) are subdominant in the macroscopic large- ./V limit, 

P^^^^-y if (IT7)2 + (1^7)2 - 1 ' ^ 
where z = x + iy. The density is constant and bounded by an ellipse, where for r = we recover the 



result of [14]. On the other hand for eq. (2^8) the macroscopic density has been calculated in 



1 / fi 2 + y 2 



where the eigenvalues are bounded by the following curve) 2 ] 

= y V + 4// 2 ) + y 2 (l + 2// 2 (2 - x 2 - 4/i 2 )) + //(x 2 - 4(1 - fi 2 )) . (2.11) 



The aim of the two models (|2.l| ) and ( |2.8[) is to describe the fluctuations of Dirac eigenvalues close 
to the origin, because of their relevance for chiral symmetry breaking from the Banks-Casher relation 
p4|. We will therefore restrict ourselves to small values of fj, since otherwise the determinant in eq. 
( |2.8| ) is completely dominated by fi. Another reason for small \i is that we wish to stay in the phase 



with broken chiral symmetry. For fi > 1 the curve (2.11) splits into two and the vanishing of the 
macroscopic spectral density |0u(0) indicates through the Banks-Casher relation that the chiral phase 
transition has taken place. If we solve the biquadratic equation ( 2.11| ) for the imaginary part y and 
expand in \x we obtain 

y 2 = /A4-* 2 ) + O(^) . (2.12) 



Expanding the density eq. (2.10) in /i < 1 we can therefore neglect y and finally arrive at 



Comparing eqs. ( |2.9D and ( |2,13| ) we identify 

4^ = (1 _ r 2 ) j (214) 

which is valid for small chemical potential \x and the parameter r close to unity meaning small non- 
Hermiticity. We have now identified all the parameters in our model with the corresponding quantities 
in QCD. In contrast to the model ( [2.8; ) by Stephanov we will not be able to study the chiral phase 
transition since our p T (z) will remain constant for all r G [0, 1). On the other hand we will be able 
to calculate all correlation functions, for finite- iV as well in two different microscopic scaling limits in 
the next sections. This has not been achieved so far for the model ( |2.8| ), where only the macroscopic 
density is known 

We finally comment on the presence of the absolute value in the weight eq. ([T^) which is chosen 



to make the partition function real. In the model (2.8) no absolute value of the Dirac determinant is 



present. Let us however emphasize that also in the real eigenvalue model eq. ( |2.3[ ) an absolute value is 
present due to the mapping from the positive to the full real line, even for zero quark flavors Nf = 0. 

2 In our conventions the Dirac eigenvalues are real for fi = 0, as it should have become clear by now. In comparison 
to ref. |8| this means that we interchange real and imaginary part, x <-> y, compared to there. 
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3 The Solution from Orthogonal Polynomials 

We now turn to the solution of our model for a finite number of eigenvalues N, where we use the 
powerful method of orthogonal polynomials . In the two following sections we then take the large- 
N limit. The orthogonal polynomials in the complex plane we are searching for are defined with 
respect to the weight function eq. ( p.2[ ), 



dzdz* w^(z) P^ a \z)P l {a \z*) = 5, 



kl 



(3.1) 



All fc-point correlation functions defined as 

iV 



{a) t 
P N [Zl, 



. N N 

) = J H dzjdz* Hw^( Zl ) \A N (z!,...,z 2 N )\ 2 . 



(3.2) 



j=k+l 1=1 

can be expressed in terms of the kernel of orthogonal polynomials 

N-i 



(3.3) 



fc=0 



Due to the invariance of the determinant the monomial powers inside the Vandermonde A at in eq. 



( [2.1| ) can be replaced by the orthogonal polynomials and we arrive at [16] 



f$ {zi, . . . , z k ) = det 

l<i,i<fc 



(3.4) 



In eq. (2.1) only even powers of zj occur inside the Vandermonde A^v and we thus only need even 
polynomials. From the known relation between Hermite and Laguerre polynomials for a = — 1/2 
we expect to obtain Laguerre polynomials of the squared argument, pj t a \z) ~ L^,(z 2 ), and we will 



find that this is indeed the case (see eq. ( 3.15Q below). 

In order to derive the orthogonality of the Laguerre polynomials in the complex plane we define 
the normalization integral 



/W(r) = / dzdz* w^ a \z) 



(3.5) 



.V-' l -lvrr(a + ^ (l-r 2 )f+!P a . 



where P a+ 3 (x) is the Legendre function. While the expression given is exact for finite and infinite iV it 

will simplify in the weak non-Hermiticity limit in section || where 1 — r 2 is rescaled with N. Following 
the idea of ref. |J we use the invariance of the integral eq. (|3.5|) under reparametrization to make the 



generating function of the Laguerre polynomials appear. Performing a change of variables z 



e'^z 



inside the integral eq. ( |3lf ) we obtain from eq. (2.2) 



/(«)(, 



dzdz*\z 



*|„|2a+l 



exp 



N 



cxp 



1-T 2 

(e^-l)r 2 (Nz 2 



z* + e 



1-T 2 



2r 



exp 



( e -*2 V _ i) T 2 



Nz*- 
2t 



(3.6) 
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where the average is taken with respect to the weight function eq. (|2.2| ) and we have divided by f^ a \r) 
in order to normalize. We now wish to obtain the following generating functional for the Laguerre 
polynomials, which is also valid for complex arguments Z: 



(1-u) 



-a-l 



exp 



u - 1 



E 

fc=0 



L a k (Z 2 ) u k 



\u\ < 1 , 



(3.7) 



and similarly for complex conjugate arguments. The split into variables u and Z 2 is not unique in eq. 
(p.6|). We will use the known results || for the Hermite polynomials He n of the model eq. ( |2.6| ) to 
read off the argument in the case a = —1/2 using 



He 2k U^z 



(-l) k 2 k k\ L k l 



Nz 7 
"27 



(3.8) 



We therefore choose Z = and 



which can be written as 



u — 1 



(e i2 ^ - l)r 2 
= 1-r 2 

r 2 (l-e i2 ^) 



(3.9) 



(3.10) 



I _ r 2gi2(/3 

It is complex and satisfies \u\ < 1. We now multiply eq. (|3.6| ) with [(1 — n)(l — u*)] _cl_1 to obtain 

'Nz* 2 



00 / /Af7 2 \ 



2r 



(3.11) 



where we have inserted the generating functional eq. (3.7). Using the fact that 

4t 4 sin 2 if 



uu 



(1 - r 2 ) 2 + 4r 2 sin 2 ip 



(3.12) 



together with the definition (|3.9| ) it is easy to see that the left hand side of eq. ( |3. 11 ) only depends 
on the combination uu* and we obtain 



[(l-u)(l-u*)] 



a— 1 



[1 - r- W] — 1 = £ ^ 



fc=0 



r(a + l + /c) (««* 
a + l)fc! r 



2A- 



(3.13) 



which converges due to t~ 2 uu* < 1 from eq. ( |3.12| ). Consequently also the right hand side of eq. 
( |3. 11 ) depends only on the combination uu* and thus the Laguerre polynomials have to be orthogonal. 
By comparing coefficients of eqs. (3.13) and ( 3.11| ) we arrive at 



7fel 2r J k \ 2r 



>kl 



r(a+l + A:) 1 
T{a+l)k\ T 2k 



(3.14) 



The orthogonal polynomials from eq. (3.1) are thus finally given as 



/ru ( i ;;;;, j ;;' .r : v./r,; 



r(a + l + /c) , f > 



"27 



(3.15) 
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The phase factor (— l) fc which is chosen arbitrarily is taken such that the relation (|3.§| ) is reproduced 
at a = —1/2. All Appoint correlation functions are therefore completely determined inserting the 
polynomials Pu a \z) into eqs. ( |3.3[ ) and (|3.4| ), We only give the expression for the kernel at finite- N, 



N-l 



v *2 



. ^ ( ^^. (22 , 1 ^g f ^ T ^ J (^) iJ (^ j . (3 , 6 ) 

Let us point out that although we have checked our results by comparing to the orthogonal Hermite 



polynomial of the model (2.6) at a = —1/2 we cannot do the same check with the correlation functions. 
The reason is that in the model ( |2.6| ) also the odd polynomials contribute to the correlation functions. 

In the following sections we will investigate the large- N limit. For that purpose usually a different 
form of the kernel eq. ( |3.3[ ) is given. It only contains the polynomials of order N and N — 1 and is 
derived from the Christoffel-Darboux identity. However, due to the r dependence of the arguments 



in eq. ( 3.15| ) this identity no longer holds in the complex planef^. The modified Christoffel-Darboux 



identity now contains polynomials of all orders, 0, . . . , N. A similar fact occurs in the two-matrix 
model (see e.g. [|40|j ). The large- N limit is no longer governed by the asymptotic polynomials alone 
and we will have to use special properties of the Laguerre polynomials to achieve our goal in the next 
sections. 

Before we continue let us explain why we did not include Dirac quark masses into our model, in 
contrast to p3| ]. In the non-chiral model mass terms of the form Ylf* \ z ~ ^ m f\ 2 were included in the 
weight function. The corresponding correlation functions were calculated by relating them to massless 



correlations, using the fact that the masses can be produced from a larger Vandermonde |41, |23| . 
Using the same technique in our model eq. ( |2.l| ) we would only be able to treat additional terms of 
the form ri/ / (- 2;2 + m 2 ){z* 2 + Tn 2 ) in the weight ( |2.2|) . While these terms might be interesting for 
other applications they are not of the right form for QCD mass terms of a complex Dirac matrix. For 
this reason we refrain from deriving such correlation functions here. 



4 The weak non-Hermit icity limit 

The large- N limit at weak non-Hermiticity was introduced in |n| and is defined as follows. We take the 
limit r — > 1 that leads to Hermitian matrices at the same time as iV —* oo such that the combination 

lim N(l - r 2 ) = a 2 (4.1) 

N— >oo 



is kept finite. In this limit the macroscopic density eq. (|2.9| ) depends on real eigenvalues only. It 
becomes the Wigner semi-circle, p(x) = — x 2 , for the Gaussian weight we have chosen, with 

support [—2, 2] on the real line only. However, the properly rescaled microscopic correlations remain 
complex and are different from those of real eigenvalues as well as from the strong non-Hermiticity 
limit where r < 1. Because of the identification we have made between r and //, eq. ( |2.14 ), eq. (4.1) 
invokes to send also /i — > as N — > oo while keeping 

ANfj, 2 = a 2 (4.2) 

fixedQ. In other words the weak non-Hermiticity parameter a 2 directly measures the chemical potential 
in the microscopic scaling limit. The rescaling of \x is similar to that of the quark masses rrif, however 

3 I wish to thank E. Kanzieper for pointing out this fact to me. 

4 In ref. 0] an attempt was made to relate /i to q by comparing the macroscopic density (2.10) to the average of 
microscopic eigenvalues and a wrong power was obtained. 
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with a different power in N. It makes sure that the eigenvalues in the Dirac determinant are not 
completely dominated by the masses or the chemical potential (see in eq. ( |2.8| )). Namely in the weak 
non-Hermitian limit at the origin we also rescale the complex eigenvalues keeping 



N($le z + i9m z) 



Nz =£ 



(4.3) 



The matrix size N corresponds to the volume on the lattice. It has been mentioned in [32] that in the 
RMM eq. ( |2.g| ) the numerical effort to obtain convergence grows exponentially with Np?. Keeping it 
fixed here according to eq. ( |4.2| ) should make a comparison to data feasible. 

Eqs. ( ^4.1|) - (^^) define our microscopic origin scaling limit in the complex plane at weak non- 
Hermiticity. The kernel eq. ( |3.3| ) and correlation functions eq. ( |3.4[ ) also have to be rescaled with 
the mean level spacing 1/N of the eigenvalues. In the definition of the correlation functions for each 
variable Zj we integrate over its real and imaginary part and we obtain one power 1/N from each 
integration. The microscopic kernel and correlation functions are therefore defined as follows in the 
weak non-Hermiticity limit: 



lim 

lim 

N—foo 



— K {a) 

N 2 K N 



1 



N2k 



(a) 

Pn 



6 8 
N' N 

N''"' N 



(4.4) 



When looking at the kernel at finite- N eq. (3.16) and the appropriate rescaling of the eigenvalues, eq. 
(|4.3|), we observe that we are missing a power in N inside the argument of the Laguerre polynomials. 
The asymptotic of the Laguerre polynomials, 



lim k- a L a k {Z 2 ) 

k^oo 



(kZ 2 )- ! lj a (2VkZ 2 ~) 



(4.5) 



suggests to introduce a variable t = jj and thus to replace the sum in eq. ( |3,16j ) by an integral, 
^2k=o ~ * N Jo dt- In doing so we write everything in terms of the variables t and £j. Under the 
integral we can thus replace everywhere the Laguerre polynomials by J-Bessel functions using eg. 
( |4.5| ) since a finite t = G (0, 1] means that k has to become large and the point t = is of measure 
zero. 

In the following we separately give the asymptotic of each factor in eq. ( |3,16 ) before stating the 
final answer. The weight functions contain already the right powers of iV in the exponential and using 



eq. (fOJ) we obtain 



lim w {a) (z)^ = iV- a -5|£| a +5 exp 



1 



a- 



(4.6) 



To calculate the weak limit of the normalization integral /^(t) in eq. ( |3.5| ) we need to know the 
asymptotic of the Legendre functions, given by 



lim P ' i 

T _»l a +2 



2 a+ 2r(a + 1) 

v^r(a + §) 



(1 



This leads to 



a _ 1 
' 2 4 



lim/( a )(r) = V2^ 2 ~N- a - 2 2 a T{a + 1) . 

T— *1 



Under the integral we first replace 



lim —, ■ — 

fc^oo T(a + 1 + k) 



N~ a t- 



(4.7) 
(4.8) 

(4.9) 



10 



Rewriting eq. ( [Q| ) as r = 1 - a 2 /N the power in r turns into an exponential, 



lim ( r 

N^oo 



2k 



3 fcln(-r 2 ) 



Finally The Laguerre polynomials are replaced as 

'Nz 2 ' 



lim L% 



(4.10) 



(4.11) 



Putting together eqs. Q4.6D - ( ^4.11 ) all powers of N cancel as they should and we obtain 



4 a) (*i.Q) = l^i^e-*^^^') ^ e -° 2 V a (^i)Ja(^). (4.12) 

When inserting the kernel into the determinant eq. (^J) for the correlators all factors apart from the 
integral can be taken out. In particular the second factor cancels due to the occurrence of the £~ a 
and £j~ a in each column and row. We arrive at the microscopic, weakly non-Hermitian correlation 



functions 

» 



n 



f = \ VV2vra 2 



e a 



det 

i<*ii<* 



dt e~ a t J a (V2tCi)J a (V2t^ 



(4.13) 



They are clearly different from the correlations of the non-chiral model eq. Q2.6| ) as calculated for 
a = in |l9| and for a = Nf massless flavors in p3| . There, the product of two J-Bessel functions is 
replaced by a single cosine, cos(\/2i(£j — £p) which corresponds to a half-integer Bessel J_i. As an 



A6 




0.4 

9m £ 



Figure 1: The microscopic density at a 2 = 0.6 for a = (left) and a = 1 (right) massless flavors. 



important check we can compare to the correlation functions of real eigenvalues [24] of the chGUE eq. 
( [2.3; ) which have been shown to be universal ^3|. We take the Hermitian limit by sending a 2 — > 
which corresponds to r = 1. The exponential prefactors become 5- functions for the imaginary parts of 
the eigenvalues, as they should. The integral inside the determinant then makes the universal Bessel 
kernel (of squared arguments) appear, reproducing the known correlation functions of p4|: 



I' dtJ a (V2t Xl )J a (V2t X2 ) = E^Cy^iW^) - X2Ja + i(V2x2)J a (V2 Xl ) 
/o V2(xf - xj) 
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To give an example for eq. ( 4.13| ) we show the microscopic density in the weak limit for two values of 
a in Fig |l|. For one quark flavor (a = 1) the level repulsion at the origin is stronger as expected. The 
plot continues symmetrically from the given first quadrant into the full complex plane by reflecting 
along the real and imaginary axis. In cutting along the positive real axis we can see the oscillations 
familiar from the correlations of real eigenvalues. Each maximum corresponds to the location of a 
single eigenvalue in the complex plane. We also show the a-dependence of the quenched microscopic 
density in Fig. ||. For increasing a and thus for increasing rescaled chemical potential /i the eigenvalues 
spread further into the complex plane. In particular the density also builds up along the imaginary 
axis. As we will see in the next section in the limit a-»oowe obtain the microscopic correlations at 
strong non-Hermiticity to be defined below. The limiting picture to be compared with is Fig. ||] below 
with the quenched density on the left. 



pi 0) (£) 




.<4 t 
ym £ 



Figure 2: The quenched microscopic density (a = 0) at a 2 = 1.6 (left) and a = 2.5 (right). 



5 The strong non-Hermiticity limit 



We now turn to the strong non-Hermiticity limit. In this limit r G [0,1) is kept fixed in the large- N 
limit and consequently also /i from eq. ( [2.14 ) is kept finite and unsealed, in contrast to eq. fl4.2j ). The 
eigenvalues are now rescaled differently, namely with the square root of N pll ^3[, 



VN (5fte z + iQm z) = VN z =£ . 



(5.1) 



This defines our microscopic origin limit at strong non-Hermiticity. In consequence we also have to 
change the rescaling of the microscopic kernel and correlation functions accordingly, 



1 



(a) ( a £ 2 * 



1 (a) / £l gfc 

\ -00 N k Pn \ Jn'"' 



(5.2) 



In the strong scaling limit eq. ( |5.1[ ) the argument of the Laguerre polynomials in eq. ( p. 16 ) can be 
rewritten in terms of the scaling variable £ alone. In particular we cannot use the asymptotic eq. 
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( |j.5| ) for the Laguerre polynomials any longer. Fortunately there is an identity for an infinite sum of 
Laguerre polynomials at hand [37], 



kl 



f^T(a+l + k) feV u k\ * J 



k=0 



+ k) 



Z 2 Zfr 2 ) 

1-T 2 



■ exp 



1-T 2 



(Zf + Z 



*2 \ 



l-r 2 



(5-3) 

which is the analog of Mehlers formula for the Hermite polynomials used in [19] to derive the strong 
correlation functions. We can thus immediately read off the large- iV limit of the kernel eq. ( [3.16 ) and 
obtain 



4%, 6) 



2 a T(a + 1) 



* \a+T 



vrr(a + |) (l-r2)f + iP a+ 3 



feci) 



x exp 



2(1 -r2) 



* 2 



la 



1-T 2 



(5.4) 



Here we have inserted the explicit expression for the normalization integral eq. (3.5). After taking 
out factors from the determinant the correlation functions again simplify considerably and our final 
result reads 



( a ) 1 1- 
Ps (&> 



n 

i=i 



2 a T(a + 1) \i 



vrr (a+§)(l 



det 

l<i,j<k 



(5-5) 

In ref. [p7o| ^ correlation functions for a charged two-dimensional plasma with fractional charges at the 
origin have been obtained, in what we call strong non-Hermiticity limit. However, since the interaction 



term induced by the Vandermonde determinant of squared arguments in our model eq. (2.1) is different 
from that of ]i~0[ ] the correlation functions obtained there also differ from ours. 

A consistency check for the correlation functions eq. ([T5|) can be obtained from the weakly non- 
Hermitian correlations eq. ( f4 . 1 3| ) by taking the limit a — ► oo there. We evaluate the integral inside 
the determinant as follows, 



/ dt e~ a 1 J a (V2tti) J a {V2ie 

Jo 



1 f°° 

I dse s J a 



a 2 Jo 



a*- 



exp 



2a 2 




(£f + ^ 



(5-6) 



When identifying a 2 = 1 — r 2 and taking out the exponential of the determinant in eq. ( |4.13 ) we 



immediately recover the strong correlations eq. (5.5). In order to match the normalizations we also 
have to expand the Legendre function for large argument, according to eq. ( ^.7|) . We have thus 
control over the full parameter range starting from real eigenvalues via weak to strong non-Hermitian 
correlation functions. 

As an example for eq. ( |5.5| ) we depict the strong microscopic density in Fig. ^ for zero and one 
massless flavors. The value for r = 0.5 corresponds to fi = \/3/4 ~ 0.43 in our units, according to 
eq. ( |2.14| ). The microscopic spectral density develops a hole at the origin and becomes flat for large 
values. This has indeed been observed in quenched lattice data for intermediate values of the chemical 
potential [p3| |. For one massless flavor (a = 1) the level repulsion at the origin is again stronger and 
the barrier gets flattened out. The picture is very much reminiscent of the microscopic density with 
massless flavors of the non-chiral model [p3|[ in the strong limit (compared to Fig. 3 there). When 
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Figure 3: The microscopic density at r = 0.5 for a = (left) and a = 1 (right) massless flavors. 



getting closer to the Hermitian limit by raising r / 1 or equivalently lowering fj, \ the repulsion 
at the origin gets much more pronounced. In Fig. || we have depicted the same densities as in Fig. 
IH but for r = \/2T/5 ~ 0.9 corresponding to /x = 0.2. Both values for in Figs. || and |] are in the 
broken phase when comparing to ref. Q. There, the phase transition occurs at a critical value fi c = 1 
as it follows from eq. ( 2.1l| ). Even if we take the extreme case of maximal non-Hermiticity which 
corresponds to r = in our model we will stay inside the phase with broken chiral symmetry. This can 
be seen when assuming that the relation ( |2.14j ) is still approximately valid, leading to /i(r = 0) = ^. 




9m £ 



Figure 4: The microscopic density at r ~ 0.9 for a = (left) and a = 1 (right) massless flavors. 



14 



6 Conclusions 



We have introduced and solved a new class of chiral RMM with complex eigenvalues, that corresponds 
to the extension of the chGUE into the complex plane. It was shown that the orthogonal polynomials 
of our model are given by the Laguerre polynomials in the complex plane. Consequently all k- 
point correlation functions could be given exactly for a finite number of eigenvalues N in terms of 
the kernel of these polynomials. We then investigated two different large-iV limits with weak and 
strong non-Hermiticity. In both cases explicit results were given for all correlation functions in terms 
of the two different limiting kernels. Furthermore, we could prove that the results for weak non- 
Hermiticity interpolate between the known correlation functions of real eigenvalues of the chGUE and 
those of strong non-Hermiticity, taking the limit of the weak non-Hermiticity parameter a — > or oo, 
respectively. 

As an application of our model we have proposed that it describes the correlations of QCD Dirac 
eigenvalues close to the origin in the presence of a small chemical potential fi. For that purpose we 
related our parameter r that measures the non-Hermiticity with the chemical potential as it occurs in 
QCD. In the identification we compared to another matrix model that incorporates the symmetries 
of the Dirac operator with \i ^ more closely. It remains to compare our predictions to lattice data 
in four-dimensional QCD yet to be generated and to see if the success of RMM in predicting real 
eigenvalues extends to the complex plane. It would be also very interesting to find other applications 
too, in particular given the interpretation of the complex eigenvalues as positions of charges distributed 
in a two-dimensional plane. 

So far we have only defined our model as an eigenvalue model and obtained its chiral interpretation 
by analogy with the corresponding real eigenvalue model, the chGUE, and its non-chiral counterpart, 
the GUE. It would be very interesting to find an explicit matrix representation of our model, in 
particular in view of the definition of Dirac mass terms and the corresponding solution. 

The structure of RMM with complex eigenvalues seems to be much richer as the number of sym- 
metry classes is much larger than those of real eigenvalue models. While not all of these model may 
be accessible to the technique of orthogonal polynomials presented here at least the chiral versions 
of the non-Hermitian or non-symmetric matrices with Dyson index (3 = 4 or /3 = 1, respectively, 
could be worth investigating. The fact that our chiral model with [3 = 2 again contains the classical 
polynomials of Laguerre type may be of help. 
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